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Abstract. We study Harper operators and the closely related discrete magnetic Lapla- 
cians (DML) on a graph with a free action of a discrete group, as defined by Sunada 



Sun|. A main result in this paper is that the spectral density function of DMLs asso- 
ciated to rational weight functions on graphs with a free action of an amenable discrete 
group, can be approximated by the average spectral density function of the DMLs on a 
regular exhaustion, with either Dirichlet or Neumann boundary conditions. This then 
gives a criterion for the existence of gaps in the spectrum of the DML, as well as other 
interesting spectral properties of such DMLs. The technique used incorporates some 
results of algebraic number theory. 



Introduction 

Both the Harper operator and the discrete magnetic Laplacian (DML) on the Cayley 
graph of Z 2 have been extremely well studied in mathematical physics, as they arise as 
the Hamiltonian for the discrete model describing the quantum mechanics of free electrons 
in the presence of a magnetic field. In particular, the DML is the Hamiltonian of the 
discrete model for the integer quantum Hall effect, cf. [Bel]. These operators can be 



easily generalized to the Cayley graph of an arbitrary discrete group. This and a further 
generalization to general graphs with a free co-compact action of a discrete group with 



finite quotient, was defined by Sunada [ Bun ] and studied in the context of noncommutative 



Bloch theory and the fractional quantum Hall effect in [pHMMj , [|CHMf , fMryf , flMM2|l 



In this paper we are primarily concerned with the approximation of the spectral density 
function of the DML by examining restrictions of the DML over a sequence of finite 
subgraphs. Consider a graph X which has a finite fundamental domain under the free 
action of an amenable group T. By F0lner's characterization of amenability (see also 
|Ad|| ) there must exist a regular exhaustion of T: a tower of finite subsets A m C A m+1 , 
U m A m = T satisfying 

(0.1) km ~~7lT = 

where c^A^ = {7 E T : ^(7, A^) < 5 and £^1(7, T \ A^) < 6} is a ^-neighborhood of the 
boundary of A&, and d\ is the word metric on T with respect to some generating set. 
Given a choice 5F of fundamental domain for the V action on X, a corresponding sequence 
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of subgraphs X m is constructed, with X m being the largest subgraph of X contained 
within U{73~| 7 6 A m }. These X m satisfy a similar property to ( p.l[ ): 

(0.2) hm #VertdsX m =Q ytf> 

m~>oo ^fcVertA m 

where c^Xm refers to the subgraph which is the <5-boundary of X m in X, that being the 
intersection of the ^-neighbourhoods of X m and X \ X m in the simplicial metric. 

The DML restricted to the space of functions over the vertices of X m is finite dimen- 
sional, and these restricted DMLs constitute a sequence of finite approximations to the 
DML itself. This leads us to the following conjecture. 

Conjecture 0.1. Let X be a graph on which there is a free group action by an amenable 
group V, with finite fundamental domain. Let {-X m } oc _ 1 be the regular exhaustion of X 
corresponding to a regular exhaustion A m ofT. Then 

(0.3) lim %^1 = F(A), VAeJR, 

m^oo N m 

where N m is the size of A m , E m (\) denotes the number of eigenvalues less than or equal 
to A of the DML restricted to X m with either Dirichlet or Neumann boundary conditions, 
and F is the spectral density function of the DML on X, which is defined using the von 
Neumann trace. 

We prove two approximation theorems which partially resolve this conjecture. The first 
states that the equality (|0.3| ) holds at all but at most a countable set of points. 



Theorem 0.2 (Rough approximation theorem for spectral density functions). In the no- 
tation above, one has 

(0.4) lim %^ = F(A) 



at every point of continuity of F. In particular, ( \0-4 ) holds for all but possibly a countable 
set of points. 

Alternative proofs for this theorem exist in mathematical physics, using for instance 
the 'Shubin formula' (see [Bel2]). We thank J. Bellissard for bringing this reference to our 



attention. 

Our main result is a refined approximation theorem, that holds only for DMLs associ- 
ated with rational weight functions. 

Theorem 0.3 (Refined approximation theorem; the rational case). In the notation above, 
suppose that the DML is associated with a rational weight function. Then one has 

(0.5) lim %^ = F(A), VAel. 

m—>oo N m 



The strategy of proof of these approximation theorems is similar to that in [ DM ] , but 



the details differ and are more involved in this case. The key lemma for the proof of 



Theorem 0.2 is the combinatorial analogue of the principle of not feeling the boundary. 
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Lemmas 2.1 and 2.2 in section 2 establish an approximation for F(\) in terms of the 
and a sequence of polynomial approximations to the characteristic 



restricted DMLs 
function X[o,A] : 

(0.6) 



(m) 



F(X) = lim lim 



+oo m^oo 



Tr 



c I P 



The proof of the refined approximation theorem |0.3| has at its core a variant of Lemma 
2.8 in [Lu|. Such a result is called 'Log Holder continuity' in the mathematical physics 
literature. By providing bounds on the growth of scaled spectral densities E m (X)/N m , 
our variant of Luck's lemma allows us to effectively interchange the limits in ( |0.6j ) and 
provide us with our result. Showing that the lemma is applicable is more involved in the 
case of the DML; we require a positive lower bound on the modulus of the product of 
non-zero eigenvalues of the restricted DMLs translated by A. Such a bound for algebraic A 
is established in section 3, using some interesting facts about algebraic number fields and 



their rings of integers. Proposition 2.5, which states that the point spectrum of the DML 
with rational weight function must be contained within the algebraic numbers, completes 
the proof. 

The conjecture is known to be true for all weight functions - not just rational ones - 
when X is the Cayley graph of Z 2 . In this case the spectral density function of the DML 
can be shown to be continuous (|Sh|, JDS|]), and the result then follows from Theorem p.2| . 
The proof that the spectral density function is continuous appears to utilize the special 
geometry of the Cayley graph of Z 2 . We do not know if the spectral density function 
is continuous for general amenable graphs, and our proof of Theorem |(0| is thus quite 
different to that of the Z 2 result. 

In the final section of this paper we examine some consequences of these approximation 
theorems. Where the weight function of the DML is rational, Corollary |4.2| provides a 
criterion for the existence of spectral gaps which is in principle computable; the existence 
of spectral gaps for general DMLs is unknown and is a central open question in the area. 
Further in this section, we consider the Fuglede-Kadison determinant of DML— A, and 
show that it is greater than zero for most A. This qualification on the spectrum of the 
DML may have applications in Mathematical Physics. In joint work with Dodziuk, Linnell 



and Schick |DLMSY], we have also worked out the analogs of the results of this paper for 
the combinatorial Laplacian on L 2 cochains on covering spaces of finite CW complexes. 



1. The Discrete magnetic Laplacian 

Consider a combinatorial graph X on which there is a free action by a group T, with 
finite fundamental domain. The edge set EdgeX is a collection of oriented edge; each 
combinatorial edge in X has two corresponding oriented edges in EdgeX, one for each 
choice of orientation. 

If e is an oriented edge, e will denote the edge with opposite orientation, and t(e) and 
o(e) will denote the terminus and origin respectively. 

It will be convenient to regard a subset E + of Edge X in which each combinatorial edge 
has exactly one oriented representative; E + corresponds to a choice of orientation for the 
graph. 
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If g is a function defined on the edge set EdgeX with values in some group H, then 
for compatibility with the orientation reversing operation on the edges, we demand that 
g(e) = g(e)~ l . The graph structure gives a coboundary operator d from the space of I? 
functions on the vertices of X to the L 2 functions on the edges, giving the chain complex 

0^C? 2) (X)±Ct 2) (X)-^0 

where 

(df)(e) = /(t(e))-/(o(e)). 

Here t(e) denotes the terminus of the edge e, and o(e) the origin. 
The adjoint d* is given by 

(d*g)(v) = £ g(e) - £ g(e) 

eeE+ eeE+ 
t(e)=v o(e)=v 

One can then define the discrete Laplacian on (X) by 

Af = d*df 

(Af)(v) = 0(v)f(v)- £ Me))- £ /(t(e)) 

e<=E+ eeE+ 
l(e)=v o(e)=v 

where 0(v) is the valence of the vertex v. 

Comparing this with the random walk operator on C9 2 JX) 

(Rf)(v)= £ /(0(e)) + £ /(t(e)), 

t(e)=u o(e)=i> 

we see that 

(Af)(v) = Q(v)f(v)-(Rf)(v). 

The random walk operator is the basis of the generalized Harper operator of Sunada, 
where the value of / at each of the neighbouring points is weighted with a complex number 
depending on the edge. 

More precisely, consider the cochain of U(l) valued functions on the graph X. Two 
functions in C 1 (X; U(l)) are said to be equivalent if they belong in the same cohomology 
class. The weight function a that we choose on the graph must be weakly I '-invariant, 
meaning that a is equivalent to its left translation by any element of T. Specifically, for 
each 7 e T there must exist a function s 7 : Vert X — ► U(l) such that 

cr( 7 e) = cr(e)s 7 (t(e))s 7 (o(e)). 
The Harper operator is then given by 

(H a f)(v) = £ <7(e)/(o(e)) + £ ^)/(t(e)) 

e€E+ e£E+ 
t(e)=v o(e)=v 

and one can define the discrete magnetic Laplacian (DML) in terms of H a , 
(1-1) (A a f)(v) = 0(v)f(v)-(H a f)(v). 
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The discrete Laplacian arises naturally from the cochain complex on X. The DML can 
also be formed this way by examining a twisted complex with coboundary operator d T 

d T : C( 2 )(X) — >C(2)(X) 



K/)(e)=r(e)/(t(e))-r( e )/(o(e)). 

where r is a weight function satisfying 

While the normal discrete Laplacian commutes with the the left T-translations, the 
DML does not. In order to discuss its spectral properties we need to determine the 
appropriate von Neumann algebra in which to examine it. 

The magnetic translation operators are the maps 

{T 1 f)(x) = s 1 { 1 - 1 x)f{T 1 x). 
The form a projective representation of V with 

(1.2) TjTy = 9( 7 , 7') r 77' V 7, 7 G T, 

where is a [/(l)-valued group 2-cocycle determined by a up to cohomology. 

The DML commutes with these operators, and the natural context for the examination 
of its spectral properties is the finite von Neumann algebra B(l 2 (X)) r,cr consisting of all 
the bounded operators that commute with the T 7 . Hereafter Trr jCr will denote the unique 
trace on this algebra, with associated von Neumann dimension dimr^. 

More precisely, let A £ B(l 2 (X)) r,a and k(vi,V2), Vi G VertX denote its kernel or 
matrix. It follows that 

(1.3) Sr y (vi)k("fVi,"fV2)s 1 (v 2 y 1 = k(vi,v 2 ) 

for all 7 G r and for all v\,v<2 S VertX. The von Neumann trace of A is given by the 
expression 

(1.4) Tr r ^(A) = ^2k(v,v), 

where 3 denotes a fundamental domain for the action of T on X. The expression is well 



defined in view of (1.3) 



To express the DML in terms of an equivariant coboundary operator, it is necessary to 
define first compatible magnetic translation operators on C^(X). 

In order to define such translations that commute with the twisted coboundary operator, 
a further condition must be imposed upon s 7 . As the sole restriction on s 7 is that ds^a = 
7*<7, one can scale it by /c 7 , constant on connected components of X. The expression 
s fl (x)s~ / (nx)s~ /fM (x) is constant on connected components of X ([ fSun| ]) as is si(x). So 
setting A: 7 (x) 2 = s 7 (x)s 7 -i (72;) and letting s 7 (x) = s 7 (x)/c 7 (x), one has ds 7 = ds'^ and 
s 7 (x) _1 = s'^_ 1 ( / yx). Hereafter it will be assumed that s-y(x) enjoys this normalization 
property. 

Recall that the weakly T-invariant weight function r satisfies r 2 = a. Then one can 
choose i 7 with t 2 = ~s^ and r(je) = r(e)i 7 (t(e))t 7 (o(e)). t 7 can be extended to edges by 



t r (e)=t r (t(e))t y (o(e)). 
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Extending the magnetic translation operators to C^(X) by 

(T 75 )(e) =t 1 ( 1 - 1 e)g{ 1 - 1 e) 

then allows them to commute with the twisted coboundary operator d T , while preserving 
the property 1,2 . 

In order to obtain the bounds required for the proofs of the approximation theorems 
in the next section, we need the following lemmas characterizing the near-diagonality of 
the DML and giving a common bound on the operator norm of the DML and its finite 
restrictions A^. 

Let d(vx,V2) and d m (vi,V2) denote the distance between vertices v% and V2 in the sim- 
plicial metric on X and X m respectively. Recall that d(v\, V2) is the length of the shortest 
sequence of edges {ei, . . . , e^} with o(ei) = v±, t(ejfc) = V2 and t(ej) = o(ej + i) for all i. 

Let D(v\,V2) = (A a 5 vl ,5 V2 ) denote the matrix coefficient of the DML, A CT , and D^ m '{v\, v 2 ) 

at> 5 Vl , 5 V2 > denote the matrix coefficient of the DML, A a 



Lemma 1.1. D(vi,V2) = whenever d(v\,V2) > 1 and D^ m \vi,V2) = whenever 
d m (vi,V2) > 1. There is also a positive constant C independent of V\,V2 such that 
\D(yx,v 2 )\ < C and \D^( Vl ,v 2 )\ < C. 

Let D k (v\,V2) = (A^5 Vl ,5 V2 ) denote the matrix coefficient of the k-th power of the 
DML, A£, and D^ k (v\, v 2 ) = ((a^A 5 Vl ,5 V2 \ denote the matrix coefficient of the 



ifc-th power of the DML, A^ ■ Then 

D k (vi,v 2 )= D(v 1 ,u 1 )D(u 1 ,u 2 ) ■ ■ ■ D(u k -i,v 2 ) 

ui ,...Ufc_i GVert X 

and 

£)M*( UljU2 )= D^\v l ,u l )D^\u l: U2)...D^ m \u k ^V2) 

ui,...u k _ 1 GVertX m 

Then the following lemma follows easily from Lemma 

Lemma 1.2. Let k be a positive integer. Then D k (v\,V2) = whenever d(v\,V2) > k 
and D^ ,k (vi,V2) = whenever d m (vi,V2) > k. There is also a positive constant C 
independent ofv\,V2 such that \D k (vi, f 2) | < C k and \D^ k (vi, v 2 )\ < C k . 

Since the DML A k commutes with magnetic translations, A k belongs to B(l 2 (X)) r,a 
and the von Neumann trace is given as in ( |1.4| ), 

(1-5) Tr r , a (A k ) = ^2D k (v,v), 

Lemma 1.3. There is a positive number K such that the operator norms of A a and of 
Act™' 1 for all m = 1, 2 . . . are smaller than K 2 . 



Proof. The proof is similar to that in Lu], Lemma 2.5 and uses Lemma 1.1 together with 
uniform local finiteness of X. More precisely we use the fact that the valency of any vertex 
in X is uniformly bounded, say < b. A fortiori the same is true (with the same constant 
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b) for X m for all m. We now estimate the £ 2 norm of A ct k for a function k = a v 5 v 
Now 



it \ v / 



so that 



E 



E 1 ^ 



<E E i^)i 2 E 

« \d(w,«)<l / \d(u,v)<l 

^E E h 2 > 



a,,' 2 



< 

« d(u,u)<l 



where we have used Lemma 1.1 and Cauchy-Schwartz inequality. In the last sum above, 
for every vertex u, \a u \ 2 appears at most b times. This proves that ||Ao-«;|| 2 < C 2 6 2 ||k|| 2 . 
Identical estimate holds (with the same proof) for which yields the lemma if we set 
K = VCb. □ 

Let {E(A) : A G [0, oo)} denote the right continuous family of spectral projections of 
the DML A CT . Since A CT commutes with magnetic translations, so does E(\) = X[o,a](^<t)i 
for any A G [0, oo). Let F : [0, oo) — > [0, oo) denote the spectral density function, where 

F{\) = Tr r , CT (E(A)). 

Observe that the kernel of d is also given by E(0), and we have the following. 

Lemma 1.4. Let T be an infinite group and X a connected graph. Then 

= F(0) = dim(/eer(A CT )). 

Proof. Observe that / G L 2 (X) is in the kernel of the DML A CT = d*d T if and only if it is 
in the kernel of the twisted coboundary operator d T . That is, 



= (d T /)(e)=r(e)/(t(e))-r(e)/(o(e)) 

for all edges e G \X\\. It follows that |/(t(e))| = |/(o(e))| for all edges e G \X\\. Since X 
is connected, it follows that |/| is a constant function on X. Since T is infinite, it follows 
that / is identically zero. □ 

Let spec(A CT ) denote the L 2 -spectrum of A a . Then one can show that 
spec(A (T ) = {A G R : F(X + e) - F(X - e) > Ve > 0} . 
One can also show that A is an eigenvalue of A CT if and only if 
(1.6) liminf{F(A + e)-F(A-e)} > 0. 

In fact if A is an eigenvalue of A CT , and P\ is the projection onto the A eigenspace, then 

P A = lim{£(A + e) - E(X-e)}. 

e— >0 
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2. Proofs of the approximation theorems 

Let E m {\) denote the number of eigenvalues \i of A^ satisfying \i < A and which are 
counted with multiplicity. 

We next make the following definitions, 

E m (X) 



(2.1) 



F m (X) 
F(X) 
E(\) 
F + (A) 
Z + (A) 



N 

lim sup F m (X) 

m^oo 

lim inf F m (X) 

lim F(X + 5) 

lim F(X + 5). 

8^+0 



The following lemma can be regarded as a combinatorial analogue of the principle of 
not feeling the boundary, cf. DM |. 

Lemma 2.1. Let T be an amenable group and let p(X) = Ylt=o ar ^ r ^ e a polynomial. 
Then, 

Zh>(p(A„)) = I™ Tr c (p(A^)) . 

Proof. First observe that if v £ VertX m is such that d(v,dX m ) > k, then 
D k (v,v) = (A^„A) = (aW^,« ( ) = D^ k (v,v). 

By O) 

i V77- 



Therefore we see that 



Tr r , CT (p(A CT ))-^Tr c (p(A( m ) 



< 



u 6 Vert Xn 



l^E^I E (\D r (v,v)\ + \D^ r (v,v)\). 

m r =0 v&VertX m 



VGVertXrr 

d(v,dX m )<d 



Using Lemma [L^, we see that there is a positive constant C such that 

d 

< 



Tr ri >(AJ)-^Tr c (p(A( m ) 



2 ^#VertgA kir , 



r=0 



The proof of the lemma is completed by taking the limit as m — > 00, as the X m form a 
regular exhaustion of X. □ 



We next recall the following abstract lemma of Liick [Lu which is proved using the 
Lebesgue dominated convergence theorem. 
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Lemma 2.2. Let p n {l*) be a sequence of polynomials such that for the characteristic 
function of the interval [0, A], X[o,A](A i ); an d an appropriate real number L, 

lim p„(/i) = X[ ,a](aO and \p n {p)\ < L 

n— >oo 1 ' 

holds for each [i € [0, ||Ao-|| 2 ]. Then 

lim Trr,a(p n (A*)) = F(X). 



Lemma [2.1| and Lemma give 

F(A) = lim TV r , CT (Pn(A CT )) = lim lim — TrcfpjA^ 

If one could interchange the two limits on the right-hand side above, the proof of Theorem 
0.3 would be complete. The following lemma is a variant of Lemma 2.8 in [Lu] and will 



be used in this paper to justify the interchange of the limits. Such a result is called 'Log 
Holder continuity' in the mathematical physics literature, and was proved for the case 
r = Z k by Craig and Simon p§. 



Lemma 2.3. Let G : V — > V be a self-adjoint linear map of the finite dimensional Hilbert 
space V . Let p(t) = det(i — G) be the characteristic polynomial of G. Then p(t) can be 
factorized as p(t) = t k q(t) where q(t) is a polynomial with q(0) ^ 0. Let K be a real 
number, K 2 > max{l,||G||} and C > be a positive constant with \q(0)\ > C > 0. Let 
E(X) be the number of eigenvalues r/ of G, counted with multiplicity, satisfying r/ < A. 
Then for < e < 1, the following estimate is satisfied. 

{£(6)} - {g(0)} < -logC | log if 2 
dime V ~ (— log e) dime V — log e 

Applying this lemma to the self adjoint operator G — A, we obtain the following estimate 
for < e < 1, 

{E(X + e» - |£(A)} < - log C(A) | log(^ 2 + A) 
dime V ~ (— log e) dime V — loge 

where C(A) is the lower bound for |<?a(0)|, where p\(t) is the characteristic polynomial of 
G — A, and p\(t) = t k ^q\(t) is its factorization such that q\(0) ^ 0. 

To use this lemma, observe that A^ can be regarded as a matrix with entries in the 
ring of integers of the number field generated by the weight function, whenever the weight 
function is rational. Therefore we have the following 

Lemma 2.4. Let a be a rational weight function, that is a n = 1 for some positive integer 
n. If X is a non-negative algebraic number, and p m> \(t) = det(i — (A^ — A)) is the 
characteristic polynomial of A ™ — X, and p m ,\(t) = t k Qm,x(t) with q m> \(Q) / 0, then 
there exist constants h and Q independent of m such that \q m: \(0)\ > (aN m ) Q a m , 
where a is the size of the fundamental domain 9\ 

Proof. See section ||. □ 



In order to prove Theorem 0.3 for all A, we will also need the following proposition 



characterizing the point spectra of the DML. 



10 VARGHESE MATHAI AND STUART YATES 



Proposition 2.5. The point spectrum of the DML A a is a subset of the union of the 
spectra of the restricted DMLs A^ . In particular when a is a rational weight function 
if there exists A 6 spec point A a , then X is algebraic. 



Proof. Suppose A is not an eigenvalue of any of the A^ , that is A 6 R \ U m spec A^ . 
Let P\ be the projection onto the eigenspace for A of A a . Let X^ = d\X m f] (X \ X m ) 
be the outer 1-boundary of X m and let P m and P^ be the projections onto the functions 
with support in X m and X^ respectively. One has 

Trr>P A = -^Ti c P m P x 

< — dim Im P m P\ as ||P m P A || < L 
1 'm 



giving 



Trr aP\ < liminf dimImP m Px. 

m->oo N m 



Note that by lemma 1.1 , P m A a differs from A^ 1 P m only on d\X m . One can write 
PA —PAP 4- B P° 

— A M P 4- B P° 

where B m : C°(A£J — > C°(X m n diX m ) encodes this difference. 

Consider / G ImP m P A , / = P m g, g an eigenfunction of A^ for A. From A a g = Xg we 
have 

P m A a g = \P m g, 

and so 

(2.2) (A^P m + B m P^)g = XP rn g. 

Consider two solutions g\ and g2 of this equation, with gi\x^ = 92 \ x^- Then B m P° n (gi — 
g 2 ) = and we have 

A i ^P m {g l -g 2 )=XP m {g l -g 2 ). 

As A spec A^ , this implies that P m (gi — 52) = 0. So the value of g on uniquely 
determines / = P m g, giving 

dimImP m P A < dimC°(A^) 

and thus 

Tr r ,aP\ < hm -L#3iX m 

= 0, 

demonstrating that A is not in the point spectrum of A a . As the restricted DMLs A^ 
are all finite operators with matrix elements belonging to the set of algebraic numbers, 
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the union of their spectra must in turn be a subset of the algebraic numbers; this result 
therefore implies that any A in the point spectrum of the DML must be algebraic. □ 

Remark. This lemma is easily generalized to any local operator A in the von Neumann 
algebra whose components A(u, v) = (A5 U ,5 V ) are algebraic, where we regard as local an 
operator A for whom A(u, v) = for all u and v with d(u, v) > k for some fixed k. 



Theorem p.2| follows immediately from part (i) of the following theorem, and Theorem 



D.3 follows immediately from part (ii) and the preceeding proposition 2.5 



Theorem 2.6. Let T be an amenable discrete group. 



(i) In the notation of one has 

F(X) = F + (X) = F + (X). 

(ii) If in addition, one assumes that the weight function a is rational, then 

F(X) = F(X) = F(X), 
whenever X is an algebraic number. That is, under these assumptions, one has, 

F(X) = lim F m (X). 

m-^oo 

Proof. Fix A > and define for n > 1 a continuous function f n : M — > R by 

< 1 + i if /x < A 

1 + I - n(/i - A) if A < fi < X + I 



/nO) 



v 1 if X + - < u 



Then clearly X[ ,a](aO < /n+i(/^) < /n0«) and /n(/") -> X[0,A](m) as n ^ cxd for all ^ G 
[0, oo). For each n, choose a polynomial p n such that X^A]^) < Pn(fJ-) < fn(^) holds for 
all /i G [0, if 2 ]. We can always find such a polynomial by a sufficiently close approximation 
of fn+i- Hence 

X[o,A]0) <PnO) < 2 

and 

lim = X[o,A]0) 



for all G [0, if ]. Recall that E m (X) denotes the number of eigenvalues /U of A^- 



(m) 



satisfying fi < X and counted with multiplicity. Note that 1 1 A^ 1 1 < K 2 by Lemma |L3 . 
^Tr c (^(AM)) = J 

Jv m iv m 



' Tr c (^(A( m ))) = ^ E R»M 



£ (p»(m)-i)+ E 

jU€[0,A] ^£(A,A+l/n] 



+ E ^o*) 

^(A+l/n,^ 2 ] 
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Hence, we see that 

(2-3) F W (A) = %^ < -L TccMAi™))). 

In addition, 

Tr c (p„(A( m ))) < %^ + -L SU pK( M ) - 1 : fx G [0, A]} E m (X) 

+ — sup{p„(/x) [A, A + 1/n]} {E m (X + 1/n) - # m (A)) 

+ — sup{p n ( / u) : // G [A + 1/n, ET 2 ]} (E m (K 2 ) - E m (X + 1/n)) 

< E m (\) | Efn{X) | (l + l/n)(£ m (A + l/n)-£ m (A)) 
iV m nA?" m iV m 

W 2 )-4(A + l/n)) 



niV m 



< E m {\ + 1/n) 1 £ m (K 2 ) 



N m n N m 

< F m (X + l/n) + - 
n 

since E m (K 2 ) = dim C°(AT m ) = aN m for a positive constant a independent of m, a being 
the size of the fundamental domain 3'. It follows that 

(2-4) Tr c (p n (A^))) < F m (A + 1/n) + -. 

iv m n 

Taking the limit inferior in ( |2.4D and the limit superior in (|2.3|) , as m — ► oo, we get that 

(2.5) F(X) < Tr Tja (p n (A a )) < F(X + 1/n) + -. 

' n 

Taking the limit as n — ► oo in (^1]) and using Lemma |2,2| , we see that 

F(A) < F(A) < F+(X). 

For all e > we have 

F(X) < F+(X) < F(X + e)< F(X + e) < F(X + e). 
Since F is right continuous, we see that 

F(A) = F + (X) =F+(X) 

proving part (i) of Theorem [2.6| . 

Next we assume that a be a rational weight function, that is a n = 1 for some positive 
integer n. In the notation of Lemma 2A, let A be a non-negative algebraic number, 
and p mi \(t) = det(i — (A^ — A)) be the characteristic polynomial of A^ — A, and 
Pm,\(t) = t k q m , x (t) with q m>x (0) + 0. Then |g m , A (0)| > (aN m )- h Q- haN ™ , h and Q being 
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constants independent of m. By Lemma |2.3| and the remarks following it, for e > one 
has, 

F m (X + e) - F m (X) ^ - log((aN m )- h Q- haN ^) log(K 2 + A) 



< ^ ^ '- + 

(-loge)aiV m -loge 

/i(log aN m + aN m log Q) a log(K 2 + A) 



That is, 

(2.6) F m (\ + e) < F m (X) . 

(loge)iV m loge 

Taking limit inferior in Q2.6| ) asm->oo yields 

F(X + e)< F(X) - - QlQg(K2 + A) . 

log e log e 

Passing to the limit as e ^ +0, we obtain F_{X) = F_ + (X). A similar argument establishes 
that F(X) = F + (X). By part (i) of Theorem |J, one has F + (A) = F(X) = F + (A). 
Therefore we see that 

F(X) = F + (X) = F(X) = F + (X) = F(X) 



which proves part (ii) of Theorem |2,6| . □ 
3. Determining a lower bound for q m! \(0) 



That we can determine the lower bound in lemma 2.4 depends upon the components 



of the matrix taking values in a ring of algebraic integers. Before providing a proof 
of the lemma, some definitions and elementary results of algebraic number theory are 
presented. 

The weight function a takes values in U(l). When it is rational, that is there is some 
integer n for which a n = 1, these values are nth roots of unity. The matrix elements in 
are integer- linear combinations of the values taken by a, and so lie in where C 
is a primitive nth root of unity. 

The natural context for examining the properties of this matrix is the algebraic number 
field. 

Definition 3.1. An extension F of Q of finite degree is an algebraic number field. The 
degree of the extension \F : Q| = dimQ F is called the degree of F. The integral closure of 
Z in F (that is, the subring of F consisting of elements a that satisfy f(a) = for some 
monic polynomial / £ Z[X]) is termed the ring of (algebraic) integers of F, and written 
F . 

For a given algebraic number A, let A = r\jb where r\ is an algebraic integer and b E 7L. 
The field in which the DML will be examined will then be Q(CiV)- The matrix elements 
of A^ are elements of Z[£] and thus reside within the ring of integers of this field. Let 
h be the degree of this extension, which will depend on the denominator n of the rational 
weight function a and this algebraic integer r\. 

The extension Q(C, n) is Q(a) for some algebraic a. The minimum polynomial of a will 
be of degree h. The other roots of this polynomial - the conjugates of a - determine the 
set {ej} of Q-preserving embeddings of Q(a) into C. 

The set of embeddings of an algebraic number field determine a norm on that field. 
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Definition 3.2. Let F be an algebraic number field of degree h, with distinct embeddings 
ei, B2, ■ ■ ■ ,&h into C. Then the norm on F is defined as the product 



h 



N(a) = Y[ei(a). 



i=l 



The norm has the following important property: 

Theorem 3.3. Let N be the norm on an algebraic number field F. Then N(a) E Q for 
all a E F , and further, N(ot) E Z for all a in the ring of integers Of- 

In particular, N(a) E Z for all a E Z[£]. 



From this one can arrive at the following well-known result, presented in [Fa| and 
repeated here. 

Lemma 3.4. Let F C C be an algebraic number field of degree h, with ring of integers 
Of and distinct embeddings e±,e2,... , into C. Let a be a non-zero element of Of 
satisfying 

\ei(a)\ < R Vt. 

Then 

\ei(a)\ > R 1 - 11 Vi. 



Proof of lemma |2.4| . We follow here an argument of Farber |Fa|| . As discussed above, 
let {ej} be the h Q-preserving embeddings of Q(C,v) = Q( a ) i n to C, and let e\ be the 
identity embedding. 

Now Q = ri(a) for some polynomials r% of degree less than h. There will be an upper 
bound on these co-efficients, as the number of distinct is hnite. As there is also a bound 
on the absolute values of the conjugates of a, there is an upper bound R: 

\ej(C)\<R Vi,j 
Mr?) | < R Vj 

Let the matrix Aj = ej(A^), applying ej component- wise. Then the matrix elements 
of Aj satisfy the same conditions as those of the DML in Lemma |1.1| (for a different 

constant C) and by the same procedure by which a bound K 2 was determined for || || , 
one can find L such that \\Aj\\ < L 2 , and thus 

|Tr(4)| ^aNrJAjW 
< aN m L 2i 

The rth symmetric polynomial of eigenvalues of Aj is therefore bounded, 

\s r {A 3 )\<^ Nm+ r r ~ 1 )-^ 

< 2 2aN ™-\L 2r (as r < aN m ) 

^ ^2aN m j2aN m _ r^2\aN m 



APPROXIMATING SPECTRAL INVARIANTS OF HARPER OPERATORS ON GRAPHS 



15 



As the matrix elements of are algebraic integers in Q(a), so are the coefficients q of 
the characteristic polynomial p m (t). One has |ej(cj)| = |sj(A,)| < {AL 2 ) aNm . 

Consider p m (t + A) = q m) \{t)t k . The coefficient of t k in p m (t + A) then is equal to 



i,a(0). 



Pm (t + A) = c aNm {t + \) aN ™ + c aNm ^(t + A)^-- 1 + • • • + c 



aNm-l 



aN m )c aNm \ t aN ™ + 



,(aN m \ ! / 



aiV m - 1\ 



l^-iA )^- 1 ^- 
aN m \ x 2 / aiV m - 1A . 1 / aA m - 2\ 







c aAr m A 



1 |C a AT m -lA + I Q Jc aA r m _ 2 A Jt 



where c a jv =1. So 



aN m —k 



q m ,x(0)= E 



i=0 



i + fe 



Q+fcA\ 



Recall that A = r\jb with rj an algebraic integer in Q(a). 6 a7Vm g m ,A(0) is then an integer 
linear combination of the q and powers of 77 which makes it an algebraic integer in Q(a). 
As 2 aNm > ( aN k m ) > { aN ^~ l ), we get the following upper bound on \ej(b aNm q m , x (0))\: 

/aN m -k . , \ 

A e n iw^ 



e# aA S m ,A(0))| 



i=0 



k 



aN m -k /• 1 1 

E it V^ws^ 



i=0 
aN m —k 



i=0 

< (8RbL 2 ) aNm aN m . 

By lemma 3.4 then, 

|e# aA S m ,A(0))| > (8^6L 2 )-^(aiV m )- /l Vj 

and in particular 

|fo aA S m ,A(0)| = |ei(6 a ^g m ,A(0))| > (8RbL 2 )~ haN - (aN m )- h , 

giving 

|?m,A(0)| > 6- aJV "(8126L 2 )-' wJV '» (aN m )- h 
> (8Rb 2 L 2 )- haNm (aN m )- h . 
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Thus for suitable choices of constants h and Q independent of m, one can write the 
inequality as 

km,A(0)| >Q' haNm (aN m y h . 

In the case where A = a/b > is rational, one has |g m ,A(0)| > (SB 2 K 2 )~ naNm (aN rn )- n 
where B is the maximum of a and 6, and K 2 is the bound on the norm of the DML as 



given by lemma 1.3 



4. The Fuglede-Kadison determinant of DML— A 

In this section, we show that the Fuglede-Kadison determinant of DML— A is positive 
for most A. We begin though by deriving some corollaries to the main theorems in the 
previous sections, using the notation there. The following corollary should be compared 
to Proposition 2.5. 

Corollary 4.1. Let spec(A CT ) denote the spectrum of the DML A a and spec(A^) denote 
the spectrum of the DML A^ . Then one has 



spec(A CT ) C (J spec(A^ 



m>l 

Proof. Let Ai, A2 be points of continuity of the spectral density function F of the DML 
A CT with Ai < A2. Then by Theorem [T^, one has 

lim (F m (A 2 )-F m (Ai)) = F(A 2 )-F(Ai) 

m^oo 

We also notice that 

F(A 2 ) > F(X 1 ) ^ spec(A CT ) n (Ai, A 2 ) + 

and 

F m {\ 2 ) > F m (X\) <^ spec(A^) n (Ai,A 2 ) / 0. 
This immediately implies the corollary. □ 

Corollary 4.2 (Spectral gap criterion). The interval (Ai,A 2 ) is in a gap in the spectrum 
of the DML A a with rational weight function a if and only if 

lim {F m {\ 2 ) - F m (\i)) = 

m— >oo 

Proof. Notice that the interval (Ax, A 2 ) is in a gap in the spectrum of the DML A a if and 
only if .F(A 2 ) = F(X\). Then by Theorem |0.3| , one has 

lim (F m (X 2 ) - F m (Ai)) = F(A 2 ) - F(X 1 ) = 0. 

m^oo 

The converse is proved similarly to the proof above. □ 

Corollary 4.3 (Spectral density estimate). Suppose that a is a rational weight function, 
and F be the spectral density function of the DML A a . Then there are positive constants 
C and 5 such that for all e £ (0,5), one has for any algebraic fi, 

F(» + e)-F(»)<-^- e . 
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Proof. This corollary follows immediately from Lemma 2.3 and Lemma 2.4, together with 



Theorem D.3L □ 



4.1. Fuglede-Kadison determinants. Here we investigate under what conditions the 
Fuglede-Kadison determinant of A a — \i is positive, for rational weight function a. Re- 
call that the Fuglede-Kadison determinant detr jCT of an operator A with spectral density 
function G is defined by 

(4.1) logdetr,^) = / log|A|dG(A) 

J0+<X<L 



for some L > ||^4||, whenever the right hand side of (4.1) is not equal to — oo, and detr-o-(-A) 



is defined to be zero if the right hand side of (HT) equals — oo, cf. [FKj. The following 



proposition improves on the estimate in Corollary 4.3 



Proposition 4.4. Let a be a rational weight function. Then for fi satisfying any of 

1. [i spec A a , 

2. n = 0, 

3. [i algebraic and \i |J m spec A^ , 

the Fuglede-Kadison determinant detr, CT (A CT — A 4 ) i> s positive. 

Proof. Condition [|. In this instance, the operator A a — /x is invertible, and positivity of 
the Fuglede-Kadison determinant follows immediately. 

Condition |2| The argument for the positivity of log detr, CT A a exactly parallels that of 

section 4 of [pM|| , using the lower bound on the modified determinant |det'| A^ = q m) o(0) 



determined in lemma 2.4 



Condition M. This case calls for a slightly modified version of the argument found in 



section 4 of [ DM ], [BFK], but using the results in sections 2 and 3 instead. 

For notational convenience, let A be the DML offset by /x, A = A a — fj,, with von 
Neumann spectral density function G(A) = F(X + fi) where F is the spectral density 
function of A a . 

Let A^ be the restricted DML similarly displaced by fi, A^ = A^ — \i. Recall that 
the normalized spectral density functions of 

are right continuous, and note that the corresponding normalized spectral density func- 
tions of A^ will be given simply by 

G m {X) = F m (A + /i). 

Observe that G m (A) are step functions. Recall from Lemma |l.3| that there exists some K 
such that IIAo-H < K 2 , and HA^H < K 2 for all m. In the following, let L be an algebraic 
number greater than K 2 + so that \\A\\ < L and ||^4( m ^|| < L. 

Denote by |det'|j4.( m ) the modified determinant of A^ m \ that is the absolute value of 
the product of all nonzero eigenvalues of A^ . Choose positive a and b such that a is less 
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than the least absolute value of any non-zero eigenvalue of , and b is greater than the 
largest absolute value of any eigenvalue of . Then 

(4.2) -L log|det'|^)= / log|A|dG m (A). 

ly m Ja<|A|<6 

Integration by parts transforms this Stieltjes integral as follows. 

(4.3) f log|A|dG m (A) = (logb)(G m (b) - G m (-b)) 

Ja<\X\<b 



a<\\\<b 

- (log a) (G m (a) - G m (-a}) - \ 

a<\X\<b 



(log a)(G m (a) - G m (-a)) - / GmW - Gm{0) dX. 



Integrating |4.1| by parts, one obtains 

(4.4) log det r ,M) = ( lo S L ) ( G ( L ) ~ G (~ L )) 

+ lim {(-loge)(G(e)-G(-e)) - / ^'^ dx}. 

e-*0+ I. Je<\X\<L A J 

Using the fact that liminf (: ^o+(— loge)(G(e) — G(— e)) > one sees that 

(4.5) logdetr,,^) > (log L)(G(L) - G{-L)) - f ° {X) ~ G{0) dX. 

J0+<\X\<L * 

We now complete the proof by estimating a lower bound of log detr j0 -(^4) from a lower 
bound on -^S- log |det'|j4( m ). Idet'l^" 1 ) is the absolute value of the product of all the non- 
zero eigenvalues of A^ m \ and hence for a rational weight function a by lemma |2.4| (and 
using the notation described there) 

|detV (m) = IwWI 

> Q- haNm {aN rn )- h 

and thence 

(4.6) log |det'|^( m ) > -ha log Q - A log(aiV m ) 

for some positive constants Q and h independent of m. 

The following estimate is proved exactly as in Lemma 2.6 of [ DM |, so we will omit its 
proof here. 

(47) f gW-°W JA< | imint / G m (X)-G m (0) dx 

J0+<\X\<L A m-*oo y +<| A |< L A 



Combining ( p~2; ) and ( f4.3| ) with the inequality ( |4.6| ) we obtain 
[ G m (X) - G m (0) dX ^ (logL)(Gm(L) _ Gm{ _ L) ) 

J0+<\\\<L A 



lim (loge)(G m (e)-G m (-e)) + halog Q + -J- log(aiV m ). 

e^0+ iV m 
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By the condition under regard, fi spec , giving 

G m (e) - G m (-e) = F m (ji + e) - F m (fi - e) = for small e, 

and so one has 
(4.8) 

f Grn(X) - G m (0) 

Jo+<\X\<L 



A 



dX < (log L)(G m (L) - G m (-L)) + halogQ + A log(aiV„ 



From (|4.5[) , (|4.7| ) and (|4.8|) , we conclude that 
(4.9) log detr>A CT > (log L) (G(L) - G(-L)) - ha log Q 

- liminf { (log L)(G m (L) - G m (-L)) + A log(aiV m ) ) 

m->oo L 7V m J 

Now h/N m \og(aN m ) — > and by part (ii) of Theorem |2.6| 

G(L) = F(L + /x) = lim F m (L + Li) 

= lim G m (L) 

m—>oc 

(and similarly for — L), hence 

logdet ri(jJ 4 > -hlogQ, 

that is, 

det r , CT (A CT -/i) > 0. 

□ 
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